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Abstract 

We characterize boundedness of iterated commutators of multiplication by a 
symbol function and tensor products of Riesz and Hilbert transforms. We obtain 
a two-sided norm estimate that shows that such operators are bounded on if 
and only if the symbol belongs to the appropriate multi-parameter BMO class. We 
extend our results to a much more intricate situation; commutators of multiplication 
by a symbol function and paraproduct-free Journe operators. We show that the 
boundedness of these commutators is also determined by the inclusion of their 
symbol function in the same multi-parameter BMO class. In this sense the tensor 
products of Riesz transforms are a representative testing class for Journe operators. 

Previous results in this direction do not apply to tensor products and only to 
Journe operators which can be reduced to Calderon-Zygmund operators. Upper 
norm estimate of Journe commutators are new even in the case of no iterations. 
Lower norm estimates for iterated commutators only existed when no tensor prod¬ 
ucts were present. In the case of one dimension, lower estimates were known for 
products of two Hilbert transforms, and without iterations. New methods using 
Journe operators are developed to obtain these lower norm estimates in the multi¬ 
parameter real variable setting. 
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1 Introduction 


As dual of the Hardy space the classical space of functions of bounded 
mean oscillation, BMO, arises naturally in many endpoint results in analysis, 
partial differential equations and probability. When entering a setting with 
several free parameters, a large variety of spaces are encountered, some of 
which lose the feature of mean oscillation itself. We are interested in charac¬ 
terizations of multi-parameter BMO spaces through boundedness of commu¬ 
tators. 

A classical result of Nehari [26] shows that a Hankel operator with anti-analytic 
symbol b mapping analytic functions into the space of anti-analytic functions 
by / •-->■ P_bf is bounded with respect to the norm if and only if the symbol 
belongs to BMO. This theorem has an equivalent formulation in terms of the 
boundedness of the commutator of the multiplication operator with symbol 
function b and the Hilbert transform [H, b] = Hb — bH. 

Ferguson-Sadosky in [14] and later Ferguson-Lacey in their groundbreaking 
paper [13] study the symbols of bounded ‘big’ and ‘little’ Hankel operators on 
the bidisk through commutators of the tensor product or of the iterated form 

[H^H 2 ,b], and [Hu[H 2 ,b]]. 

Here b = b{xi,X 2 ) and the are the Hilbert transforms acting in the vari¬ 
able. A full characterization of different two-parameter BMO spaces, Cotlar- 
Sadosky’s little BMO and Chang-Fefferman’s product BMO space, is given 
through these commutators. 

Through the use of completely different real variable methods, in [6] Coifman- 
Rochberg-Weiss extended Nehari’s one-parameter theory to real analysis in 
the sense that the Hilbert transform was replaced by Riesz transforms. These 
one-parameter results in [6] were treated in the multi-parameter setting in 
Lacey-Petermichl-Pipher-Wick [18]. Both the upper and lower estimate have 
proofs very different from those in one parameter. In addition, in both cases 
it is observed that the Riesz transforms are a representative testing class in 
the sense that BMO also ensures boundedness for (iterated) commutators 
with more general Calderon-Zygmund operators, a result now known in full 
generality due to Dalenc-Ou [8]. Notably the Riesz commutator has found 
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striking applications to compensated compactness and div-curl lemmas, [3], 

120 ], 


Our extension to the multi-parameter setting is two-fold. On the one hand we 
replace the Calderon-Zygmund operators by Journe operators J* and on the 
other hand we also iterate the commutator: 


We prove the remarkable fact that a multi-parameter BMO class still en¬ 
sures boundedness in this situation and that the collection of tensor products 
of Riesz transforms remains the representative testing class. The BMO class 
encountered is a mix of little BMO and product BMO that we call a little 
product BMO. Its precise form depends upon the distribution of variables in 
the commutator. Our result is new even when no iterations are present: in this 
case, lower estimates were only known in the case of the double Hilbert trans¬ 
form [14]. The sufficiency of the little BMO class for boundedness of Journe 
commutators had never been observed. 


It is a general fact that two-sided commutator estimates have an equivalent 
formulation in terms of weak factorization. We hnd the pre-duals of our little 
product BMO spaces and prove a corresponding weak factorization result. 


Necessity of the little product BMO condition is shown through a lower es¬ 
timate on the commutator. There is a sharp contrast when tensor products 
of Riesz transforms are considered instead of multiple Hilbert transforms and 
when iterations are present. 


In the Hilbert transform case, Toeplitz operators with operator symbol arise 
naturally. Using Riesz transforms in as a replacement, there is an absence 
of analytic structure and tools relying on analytic projection or orthogonal 
spaces are not readily available. We overcome part of this difficulty through 
the use of Calderon-Zygmund operators whose Fourier multiplier symbols are 
adapted to cones. This idea is inspired by [18]. Such operators are also men¬ 
tioned in [31]. A class of operators of this type classihes little product BMO 
through two-sided commutator estimates, but it does not allow the passage to 
a classihcation through iterated commutators with tensor products of Riesz 
transforms. In a second step, we hnd it necessary to consider upper and lower 
commutator estimates using a well-chosen family of Journe operators that are 
not of tensor product type. Through geometric considerations and an aver¬ 
aging procedure of zonal harmonics on products of spheres, we construct the 
multiplier of a special Journe operator that preserves lower commutator esti¬ 
mates and resembles the multiple Hilbert transform: it has large plateaus of 
constant values and is a polynomial in multiple Riesz transforms. We expect 
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that this construction allows other applications. 

There is an increase in difficulty when the dimension is greater than two, due 
to the simpler structure of the rotation group on In higher dimension, 
there is a rise in difficulty when tensor products involve more than two Riesz 
transforms. 

The actual passage to the Riesz transforms requires a stability estimate in 
commutator norms for certain multi-parameter singular integrals in terms of 
the mixed BMO class. In this context, we prove a qualitative upper estimate 
for iterated commutators using paraproduct free Journe operators. We make 
use of recent versions of T(l) theorems in this setting. These recent advances 
are different from the corresponding theorem of Journe [16]. The results we 
allude to have the additional feature of providing a convenient representation 
formula for bi-parameter in [22] and even multi-parameter in [28] Calderon- 
Zygmund operators by dyadic shifts. 


2 Aspects of Multi-Parameter Theory 


This section contains some review on Hardy spaces in several parameters as 
well as some new dehnitions and lemmas relevant to us. 


2.1 Chang-Fefferman BMO 


We describe the elements of product Hardy space theory, as developed by 
Chang and Fefferman as well as Journe. By this we mean the Hardy spaces as¬ 
sociated with domains like the poly-disk or := for d = {di,..., dt). 

While doing so, we typically do not distinguish whether we are working on 
or T'^. In higher dimensions, the Hilbert transform is usually replaced by the 
collection of Riesz transforms. 

The (real) one-parameter Hardy space denotes the class of functions 

with the norm 

tiiviii 

j=0 

where Rj denotes the Riesz transform or the Hilbert transform if the di¬ 
mension is one. Here and below we adopt the convention that Ro, the 0*'^ Riesz 
transform, is the identity. This space is invariant under the one-parameter fam¬ 
ily of isotropic dilations, while the product Hardy space H^^(W^) is invariant 
under dilations of each coordinate separately. That is, it is invariant under a t 
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parameter family of dilations, hence the terminology ‘mnlti-parameter’ theory. 
One way to dehne a norm on is 

*=1 

Rjj', is the Riesz transform in the direction of the variable, and the 0*^ 
Riesz transform is the identity operator. 

The dnal of the real Hardy space is BMO(M^), the t-fold prodnct 

BMO space. It is a theorem of S.-Y. Chang and R. Fefferman [4], [5] that this 
space has a characterization in terms of a prodnct Carleson measnre. 


Dehne 


||^||BMO(R‘i) •— snp 

U^Rd. 


RciU essig^ 


( 1 ) 


Here the snpremnm is taken over all open snbsets U c with hnite measnre, 
and we nse a wavelet basis adapted to rectangles R = Qi x ■ ■ ■ x where 
each Qi is a cnbe. The snperscript e rehects the fact that mnltiple wavelets are 
associated to any dyadic cnbe, see [18] for details. The fact that the snpremnm 
admits all open sets of hnite measnre cannot be omitted, as Carleson’s example 
shows [2]. This fact is responsible for some of the difhcnlties enconntered when 
working with this space. 


Theorem 1 (Chang, Fefferman) We have the equivalence of norms 


~ \\b\\BMO{R'i)- 

That is, BMO(R‘^) is the dual to 

This BMO norm is invariant nnder a t-parameter family of dilations. Here the 
dilations are isotropic in each parameter separately. See also [10] and [12]. 


2.2 Little BMO 


Following [7] and [14], we recall some facts abont the space little BMO, often 
written as ‘bmo’, and its predual. A locally integrable function b : = 

X ... X Mfo ^ C is in bmo if and only if 

ll^llbmo = sup \b{x) - 6 q| < 00 

Here the Qk are dfc-dimensional cubes and bq denotes the average of b over Q. 
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It is easy to see that this space consists of all functions that are uniformly 
in BMO in each variable separately. Let Xf, = (xi,...•, , x^). 

Then b{xf,) is a function in x^ only with the other variables hxed. Its BMO 
norm in x.^ is 


||&(a^j))||BMO = sup IQ^I M \b{x) - b{xy)Qjdxy 

Qv JQv 

and the little BMO norm becomes 

||fc||bmo = max{sup ||fe(a3c)||BMo}- 

On the bi-disk, this becomes 

||6||bmo = niax{sup ||6(xi, OIIbmojSup ||6(-, X 2 )||bmo}, 

Xi X2 

the space discussed in [14], Here, the pre-dual is the space H^(T) (x) L^(T) -I- 
L^(T) ( 8 ) iL^(T). All other cases are an obvious generalization, at the cost of 
notational inconvenience. 


2.3 Little product BMO 


In this section we dehne a BMO space which is in between little BMO and 
product BMO. As mentioned in the introduction, we aim at characterizing 
BMO spaces consisting for example of those functions 6 (xi,X 2 ,X 3 ) such that 
6 (xi,-,-) and 6 (-,-,X 3 ) are uniformly in product BMO in the remaining two 
variables. 

Definition 1 Let b : ^ C with d = (di, • • • ,dt). Take a partition X = 

{h : 1 < s < /} of {1,2, ...,t} so that = {1,2, We say that 

b E if for any choices v = {ys),Vs e Is, b is uniformly in product 

BMO in the variables indexed by Vg- We call a BMO space of this type a ‘little 
product BMO’. If for any x = (xi,...,Xt) e we define x^ by removing 
those variables indexed by Vg, the little product BMO norm becomes 

\\b\\BMOx = niax{sup \\b{xi,)\\BMo} 

where the BMO norm is product BMO in the variables indexed by Vg • 

For example, when d = (1,1,1) = 1, when t = 3 and I = 2 with A = (13) and 
A = (2), writing X = (13)(2) the space BMO(i 3 )( 2 )(T^) arises, which consists 
of those functions that are uniformly in product BMO in the variables (1,2) 
and (3, 2) respectively, as described above. Moreover, as degenerate cases, it 
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is easy to see that BMO(i 2 .,.t) and BM0(i)(2)...(t) are exactly little BMO and 
product BMO respectively, the spaces we are familiar with. 

Little product BMO spaces on can be dehned in the same way. Now we hud 
the predual of BMO(i 3 )( 2 ), which is a good model for other cases. We choose 
the order of variables most convenient for us. 

Theorem 2 The pre-dual of the spaee BMO(i 2 ,)( 2 ){T^) is equal to the spaee 

) ® i' (T) + i' (T) ® <(T<W>) 

:-{/ + <,:/£ ® L'(T) and g e L‘(T) ® 


Proof. The space 

(x)Li(T) = {/ 6 L1(T3) ; HJ,H,f,H,H,f e 

equipped with the norm ||/|| = ||/||i + ||iLi/||i + ||iL 2 /||i + \\HiH 2 f\\i is a Banach 
space. Let = L^(T^) x L^(T^) x L^(T^) x L^(T^) equipped with the norm 

||(/l, /2, /s, / 4 )||wi = ll/llll + II/2II1 + ll/slll + tfilh- 

Then we see that (x)L^(T) is isomorphically isometric to the closed 

subspace 

V = {if, H,{f),H2if),H,H2if)) : / e (x) L\T)} 

of WK Now, the dual of is equal to LL® = L^(T^) x L^(T^) x L^(T^) x 
L®(T^) equipped with the norm ||(fi'i,(y' 2 ,fl' 3 ,fi' 4 )||oo = niax{||5fj||oo : 1 ^ i ^ 4} 
so the dual space of V is equal to the quotient of VL® by the annihilator U 
of the subspace V in W®. But, using the fact that the Hilbert transforms are 
self-adjoint up to a sign change, we see that 

U = {{91,92,93,94) ■ 9i + Hig2 H293 4 - H1H294 = 0 } 


and so: 

^ W^/U ^ lm(0) 

where 

^{91,92,93,94) = 9i + Hi 92 + H29Z + H1H294 
since U = ker(6*). But 

lm(0) = L®(T3) + //i(L®(T3)) + H 2 {L^{T^)) + Hi{H 2 {L^{T^))) 

is equal to the functions that are uniformly in product BMO in variables 1 
and 2. 

Using the same reasoning we see that the dual of L^(T) (x) is equal 

to L®(T3) + H 2 {L^{T^)) + //3(L®(T3)) -f i/2i/3(L®(T3)), which is equal to 
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the space of functions that are uniformly in product BMO in variables 2 and 

3. 

Now, we consider the sum’ of the spaces (x)L^(T) and L^(T) (x) 

that is 

M(am - {(/.9) : / e U(T)®Wi.(T(w))) 

equipped with the norm 

IK/) S')! = ll/llH^^(T(Ci))(g)Li(T) + ll 5 'llLi(T)(g)H^^(T(i>i))- 

We see that, if (j) : M(i 3 )( 2 ) ^ L^((T^) is dehned by (j){f,g) = f + g, then the 
image of 0 is isometrically isomorphic to the quotient of M(i 3 )( 2 ) by the space 

N = {{f,g)^ iVf(i3)(2) : / + ^ = 0} 

= {(/, -/):/£ ® i‘(T) n L‘(T) ® ffi.(T<w))), 

Now, recall that the dual of the quotient M/N is equal to the annihilator of 
N. It is easy to see that the annihilator of N is equal to the set of ordered 
pairs (0, 0) with 0 in the intersection of the duals of the two spaces. Thus the 
dual of the image of 9 is equal to BMO(i 3 )( 2 ). The norm of an element in the 
predual is equal to its norm as an element of the double dual which is easily 
computed. QED 

Following this example, the reader may easily find the correct formulation 
for the predual of other little product BMO spaces as well those in several 
variables, replacing the Hilbert transform by all choices of Riesz transforms. 
For instance, one can prove that the predual of the space BMO(i 3 )( 2 )(M'’*) is 
equal to (x) (x) . 


3 The Hilbert transform case 


In this section, we characterize the boundedness of commutators of the form 
[H 2 , [HsHi,b]] as operators on L^(T^). In the case of the Hilbert transform, 
this case is representative of the general case and provides a starting point 
that is easier to read because of the simplicity of the expression of products 
and sums of projection onto orthogonal subspaces. Its general form can be 
found at the beginning of Section 4. 

Now let b E L^(T”) and let P and Q denote orthogonal projections onto sub¬ 
spaces of L^(T”). We shall describe relationships between functions in the 
little product BMOs and several types of projection-multiplication operators. 


These will be Hilbert transform-type operators of the form P — P-*-; and it¬ 
erated Hankel or Toeplitz type operators of the form Q^hQ (Hankel), PbP 
(Toeplitz),PQ^bQP (mixed), where b means the (not a priori bonnded) mnl- 
tiplication operator on 

We shall nse the following simple observation concerning Hilbert transform 
type operators again and again: 

Remark 1 If H = P — P^ and T : L^(T”) —> L^(T’^) is a linear operator 
then 

[H, T] = 2PTP^ - 2P^TP 

and H is bounded if and only if PTP^ and P^TP are. 

Proof. 

[P - P^)T - T{P - P^) = (P- P^)T{P + P^) - (P + P^)r(P - P^) 

= 2PTP^ - 2P^TP. 


QED 


We state the main resnlt of this section. 

Theorem 3 Let b e P^(T^). Then the following are equivalent with linear 
dependence on the respective norms 

(1) b E PM(9(i3)(2) 

(2) The commutators [H 2 , [Hi, 6 ]] and [H 2 , [P 3 , 6 ]] are bounded on P^(T^) 

(3) The commutator [H 2 , [PsPi,^]] is hounded on L^(T^). 

Corollary 1 We have the following two-sided estimate 

||^||bMO(i 3)(2) ^ ||[^2, [hf3-ffl,&]]||L2(T3)^L2(T3) ^ II ^l| BM0(13)(2) • 

It will be nsefnl to denote by Q 13 orthogonal projection on the snbspace of 
fnnctions which are either analytic or anti-analytic in the hrst and third vari¬ 
ables; Qi 3 = P 1 P 3 + P^Pf-. Then the projection onto the orthogonal of 
this snbspace is dehned by = P^Ps + PiP^- We reformnlate properties 
(2) and (3) in the statement of Theorem 3 in terms of Hankel Toeplitz type 
operators. 

Lemma 1 We have the following algebraic facts on commutators and projec¬ 
tion operators. 

(1) The commutators [H 2 , [^^ 1 , 6 ]] and [H 2 , [H^jb]] are bounded on L^(T^) if 
and only if the operators PiP 2 bPf-P 2 , Pf-P 2 bPiPf-, PiP 2 bPf-P 2 , Pf-P 2 bPiP 2 
with i e {1,3} are bounded on L^(T^). 
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(2) The commutator [if 2 , is bounded on L^(T^) if and only if all 

four operators P 2 Qi 3 bQf 2 Pf-, -P2^Qi3^Qi^3^2 

are bounded on L^(T^). 

Proof. Using Remark 1 it is easy to see that 

[if2, [Hi, b]] = 4{{P2PlbP,^Pi- - - (P2^PlfcPl^P2 - P2 ^Pi^6PiP2)) 

and that the corresponding equation for [P 2 , [P 3 , b]] is also true. This, along 
with the observation that the ranges of all arising summands are mutually 
orthogonal, gives assertion (1). To prove (2) we just notice that P 1 P 3 = 
Qi 3 ~ Qn is a Hilbert transform type operator which permits us to repeat the 
above argument replacing Pi by Q 13 . QED 


The following lemma will allow us to insert an additional Hilbert transform 
into the commutator without reducing the norm. 

Lemma 2 ||P3Pf*-P2-*-6PiP2P3||2,2^i2 = ||Pf*“P2-*-6PiP211^2^2.2• 

Proof. 

The inequality ^ is trivial, since P 3 is a projection which commutes with Pf- 
and P 2 '*'. To see notice that PsPf-Pf-bPiP2P3 is a Toeplitz operator with 
symbol Pi^P2^6PiP2. So ||P3Pi^P2^6PiP2P3|| = sup,,,\\P,^P,H{; ■,x,)PiP2\\. The 
latter is just ||Pf*-P 2 '*' 6 PiP 2 ||. For convenience we include a sketch of the facts 
about Toeplitz operators we use. Let H 3 be the operator of multiplication by 
Z3, W3{f) = Z3f, acting on L^(T^). If we dehne B = Pf‘-P2-^6PiP2 as well as 

= H^3*"(P3Pi^P2^6PiP2P3)H^3” and = W^{P,^P,^P^HPiP 2 P,^)Wr 

as operators acting on L^(T^) then the sequences An and Cn converge to B in 
the strong operator topology: it is easy to see that H 3 , IT 3 *; and P 3 commute 
with Pi,P 2 ,Pi^ and P 2 "‘“. The multiplier b satishes the equation Wf'^bW^ = b 
and = Id. So we see that 

An = P^P^{Wf^P3W^)bPiP2{Wf'^P3W^). 

But ii fs then, since is a unitary operator: 

||IH3-P3lU3-(/)-/|| = ||P3iy3-(/)_PP3-(/)|| = \fp^-I){W-){f)\\ 0 (n - »), 

as tail of a convergent Fourier series. This means that VF 3 *"'P 3 H^j^ converges 
to the identity in the strong operator topology. Thus, for each / e L^(T^) we 
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have \\{An - B){f)\\ 0. So 

\\P^P^bPiP 2 \\ ^ sup \\W*^{P3P^P^bPiP2P3)W^ 

nsN 

^ \\P3P,^P,HPiP2P3l 


QED 


Now, we are ready to proceed with the proof of the main theorem of this 
section. 

Proof, (of Theorem 3) We show ( 1 ) ( 2 ) and ( 2 ) ( 3 ). 

( 1 ) ^ ( 2 ). Consider / = f{xi,X2) and g = gix^). Then [H2,[Hi,b]]{fg) = 
g- [H2,[Hi,b]]{f). So \\[H2,[Hi,b]]{fg)\\l2^r^3) = ||i^^||i 2 (T) where F{x2) = 
\\[H2, [Pfi, b]]{f)\\L2(j2y The map g Fg has L^(T) operator norm ||T||oo. Now 
change the roles of xi and x^. The Ferguson-Lacey equivalences || [H 2 , [Pfi, 6]] || ~ 
II^IIbmo give the desired result. 

( 2 ) ( 3 ). Boundedness of the commutators [H2, [Hi, 6 ]] and [H2, [H3, b]] im¬ 
plies the boundedness of the mixed commutator [H2, [H1H3, 6 ]] by the identity 
[H2, [H1H3, b]] = Hi[H 2, [H3, b]] + [H2, [Hi, b]]H3. 

( 3 ) [ 2 ). This part relies on Lemma 2. We wish to conclude from the bound¬ 

edness of [H2, [H3H1, b]] the boundedness of [H2, [Hi, b]] and [H2, [H3, b]]. To 
see boundedness of [H2, [Hi, 6 ]], let us look at one of the Hankels from Lemma 
1. Lemma 2 shows that Pf-Pf-bP 2 Pi is bounded if and only if the operator 
P^Pf Pf-bPiP2P3 is. And the latter is an operator found in the list from part 
( 2 ) of Lemma 1. The analogous reasoning shows that all eight Hankels in 1 
are bounded and so ( 2 ) is proved. QFD 


4 Real variables: lower bounds 


In this section, we are again in with d = [di,... ,dt) and a partition 
X = (/s)isssssi of {1,..., f}. It is our aim to prove the following characterization 
theorem of the space BMOi(M'^). 

Theorem 4 The following are equivalent with linear dependence of the re¬ 
spective norms. 

(1) he BMOx{^^) 
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(2) All commutators of the form ■ ■ ■, • • •] hounded in 

where kg e R and Rksjks one-parameter Riesz transform in 

direction jk^. 

(3) All commutators of the form are bounded in 

where = (jk) kei^, 1 ^ jk ^ dk and the operators Rgjis) are a 
tensor product of Riesz transforms Rgjis) = (Sikeh ^k,jk- 

Such two-sided estimates also hold in for 1 < p < oo. Remarks will be 
made in section 7. From the inductive nature of our arguments, it will also 
be apparent that the characterization holds when we consider intermediate 
cases, meaning commutators with any hxed number of Riesz transforms in 
each iterate. Below we state onr most general two-sided estimate throngh 
Riesz transforms. 

Theorem 5 Let 1 < p < oo. Under the same assumptions as Corollary 2 and 
for any fixed n = (n^) where 1 < ^ |/s|, we have the two-sided estimate 

\\H BMOxh^’^) ~ ■ ■ ■ ’ b] ■ ■ .]||i;,p(Rd)Q < \\b\\ bmOx{R<^) 

3 

where = {jk)kehj 0 < ja, < dk and for each s, there are Ug non-zero 
choices. A Riesz transform in direction 0 is understood as the identity. 

For p = 2 and n = 1 this is the eqnivalence (1) (2) and for n = 

(|/i|,..., |/;|) it is the eqnivalence (1) (3) from Theorem 4. 

Our main focus is of course on a two-sided estimate when n = (|/i|,..., |/;|) 
when the tensor product is a paraproduct-free Journe operator: 

Corollary 2 Let j = {ji,... ,jt) with 1 ^ jk ^ dk and let for each 1 ^ s ^ I, 
jis) ^ {jk)keis be associated a tensor product of Riesz transforms Rgj(s) = 
(S)kGis here the Rkjk are jjf Riesz transforms acting on functions defined 

on the variable. We have the two-sided estimate 

II^IIbMOi(R‘*) ~ sup ||[i2ij(i), . . . , [Rt j(t), b] . . .]||£p(Rd)Q < ||&||BMOi(R‘i)- 

3 


The statements above also serve as the statement of the general case for prod¬ 
ucts of Hilbert transforms. In fact, when any dk = 1 just replace the Riesz 
transforms by the Hilbert transform in that variable. In this section, we con¬ 
sider the case dk ^ 2 for 1 ^ k ^ t and thus iterated commutators with tensor 
prodncts of Riesz transforms only. The special case when dk = 1 for some k is 
easier bnt reqnires extra care for notation, which is why we omit it here. 

The proof in the Hilbert transform case relied heavily on analytic projections 
and orthogonal spaces, a featnre that we do not have when working with 
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Riesz transforms. We are going to simnlate the one-dimensional case by a two- 
step passage via intermediary Calderon-Zygmnnd operators whose mnltiplier 
symbols are adapted to cones. 

In dimension d ^ 2, a cone C c with cnbic base is given by the data ((^, Q) 
where ^ e is the direction of the cone and the cnbe Q cz centered at 
the origin is its apertnre. The cone consists of all vectors 6 that take the form 
where 9^ = {9,^} and 9j_ e 9^Q. By AC we mean the dilated cone 
with data (^, XQ). 

A cone D with ball base has data for 0 < r < n/2 and ^ e and 

consists of the vectors {?] e : d{^,r]/\\r]\\) ^ r} where d is the geodesic 
distance (with distance of antipodal points being tt.) 

Given any cone C or D, we consider its Fonrier projection operator dehned via 
Pcf = Xcf- When the apertnres are cnbes, snch operators are combinations 
of Fonrier projections onto half spaces and as snch admit nniform bonnds. 
Among others, this fact made cnbic cones necessary in the considerations in 
[18] and [9] that we are going to need. For fnrther technical reasons in the 
proof these operators are not qnite good enongh, mainly becanse they are not 
of Calderon-Zygmnnd type. For a given cone C, consider a Calderon-Zygmnnd 
operator Tc witl^ kernel Kq whose Fonrier symbol Kc e C"® and satishes the 
estimate xc ^ Kc ^ Xii+T)c- This is accomplished by mollifying the symbol 
Xc of the cone projection associated to cone C on 8“^“^ and then extending 
radially. We nse the same dehnition for Tp. 

Given a collection of cones C = (Ck) we denote by Tc, Pc the corresponding 
tensor prodnct operators. 

In [18] it has been proved that Calderon-Zygmnnd operators adapted to certain 
cones of cnbic apertnre classify prodnct BMC via commntators. As part of the 
argnment, it was observed that test fnnctions with opposing Fonrier snpports 
made the commutator large. In [9] a rehnement was proven, that will be helpful 
to us. We prefer to work with cones with round base. Lower bounds for such 
commutators can be deduced from the assertion of the main theorem in [9], but 
we need to preserve the information on the Fourier support of the test function 
in order to succeed with our argument. Information on this test function is 
instrumental to our argument: it reduces the terms arising in the commutator 
to those resembling Hankel operators. We have the following lemma, very 
similar to that in [18] and [9], the only difference being that the cones are 
based on balls instead of cubes. 

Lemma 3 For every parameter 1 ^ k ^ t there exist a finite set of directions 
Tfc e and an aperture 0 < < 7r/2 so that, for every symbol b belonging 

to product BMO, there exist cones = D{^k,rk) with ^ o^s well as a 
normalised test function f = (S>fc=i fk whose components have Fourier support 
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in the opposing cones D{—^k,rk) such that 




The stress is on the fact that the collection is hnite, somewhat specihc and 
serves all admissible product BMO functions. 

Proof. The lemma in [9] supplies us with the sets of directions as well as 
cones of cubic aperture Qk and a test function / supported in the opposing 
cones. Now choose the aperture large enough so that (1 + T)C{fk,Qk) c 
Then we have the commutator estimate 

In fact, both commutators with cones C and D are bounded and reduce to 
||^n >(^/)||2 or \\Tc{bf )\\2 respectively thanks to the opposing Fourier support 
of /. Observe that Tdbf) = TDiTdbf)) = TdTnibf)). With ||Tel 2^2 < 1, 
we see that ||T£,( 6/)||2 ^ ||rc'(b/)|| 2 - QED 

Using this a priori lower estimate, we are going to prove the lemma below. 

Lemma 4 Let us suppose we are in with d = {di, ■ ■ ■ ,t) and a partition 
X = For every 1 < A: < f there exists a finite set of directions c 

and an aperture rk so that the following hold for all b e ■ 

(1) For every 1 ^ s < / there exists a coordinate Vg e Is and a direction d ^ 
T.us and so that with the choice of cone Dy^ = D{f^y^,ryfi and arbitrary Dk 
for coordinates k e Is\{vs} and if Dg denotes their tensor product, then 
we have 

\m ,r>i, • • ■, {Ti • • •1112^2 ^ ||^||nMOi(R‘*)! 

(2) The test function f = (X)fc=i fk which gives us a large If norm in (1) has 
Fourier supports of the fk contained in D{—f^k,rk) when k = Vg and in Dk 
otherwise. 

Before we can begin with the proof of Lemma 4, we will need a real variable 
version of the facts on Toeplitz operators used earlier. 

Lemma 5 Let Dk for 1 ^ k ^ t denote any cones with respect to the vari¬ 
able. Let denote the adapted Calderon-Zygmund operators. Let K be any 
proper subset of {k : 1 ^ k ^ t}, let Dk = ^keK and To^ the associated 
tensor product of Calderon-Zygmund operators. Let be a tensor product 
of projection operators on cones D{f^kCk) or opposing cones D{—^kCk). Let 
j f K. Then 


= lFDKbPDK\L‘2{Wl)‘0- 
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Proof. 


We will establish this by composing some unilateral shift operators and study¬ 
ing their Fourier transform in the j variable. Let denote the direction of the 
cone Dj, for any I dehne the shift operator 



Si is a translation operator on the Fourier side along the direction of the 
cone Dj. It is not hard to observe that Sf = S-i. Now dehne 



We will prove that as I ^ +co, Ai ^ B in the strong operator topology. As in 
the argument in Lemma 2, this together with the fact that Si is an isometry 
will complete the proof. To see the convergence, let’s hrst remember that Si 
only acts on the j variable, and one always has the identities 


SiS-i = Id and S.ibSi = b. 


This implies 


Ai = ToAS-iTD,Si){S.ibSi)Pf,^{S.iPn,Si) 
= TnAS-iTD,Si)bPDAS-iPD,Si). 


We claim that both S-iTd.Si and S-iPdjSi converge to the identity operator 
in the strong operator topology, which then implies that Ai B as I ^ cc. We 
will only prove S-iTo -Si —>■ Id as the second limit is almost identical. Observe 
that WS-iTd.Si/ — f\\ = \\{Td. — I)Sif\\. Given any function f and any hxed 
large I S 0. Consider the / with frequencies supported in x ... x [Dj — 
l^j) X ... X 'Mf*. In this case, Sif has Fourier support in x ... x Dj x ... x 
where the symbol of T^. equals 1. Thus, for such /, we have S_iTr).Sif = f. 
The sets x ... x (Dj — l^j) x ... x R* exhaust the frequency space. With 
WTdj —/|| 2-»2 ^ 1 the operators S-iTd^Si converge to the Identity in the strong 
operator topology, and the lemma is proved. Observe that the aperture of the 
cone Dj is not relevant to the proof. QED 

We proceed with the proof of the lower estimate for cone transforms. 

Proof, (of Lemma 4) For a given symbol b e BMOj, there exist for ah 1 ^ s < / 
coordinates v = (vs),Vs ^ h and a choice of variables not indexed by Vg, 
so that up to an arbitrarily small error 


II^IIbmOx = )||bmo, 




By Lemma 3, there exist cones = D(f^y^,ry^) with directions e and 
a normalised test function fn in variables Vg with opposing Fourier support 
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such that we have the lower estimate 

where x ... x . 


We now consider the commutator with the same cones but with full symbol 
b = 6 (-,..., Due to the lack of action on the variables not indexed by Vg, in 
the commutator, we have 


[Tvi,d^^ , • • •, [Tvi,d^^ ,&]•••] ifug) - g ■ [Tvi,d^^ 


■ ■ , [Tvi,D^^,b] ■ ■ 


for g that only depends upon variables not indexed by Vg- Again using that 
multiplication operators in have norms equal to the L°° norm of their 
symbol, for the ‘worst’ L^-normalised g we have 


\\[Tvi,d^-^, ■ ■ ■, [Tvl,D^^,b] .. .](/rf5')l|L2(Rd) 

= sup ||[T^i,d„^, ..., 

X{) 

^ II 5 [^vi,D^i,b(xP)] . . .]ifH)\\L^{Rdv) 

~ )l|BMO(„j).,.(„,)(R‘i-u) = ||&||bMOi(R‘*)- 

Note that the test function g can be chosen with well distributed Fourier 
transform. Take any cones in the variables not indexed by Vg and let D denote 
the tensor product of their projections, /r = Pog- Notice that 


{Tvi,d^^ \Tvl,D■u^ ,&]•••] {IhIt) II ^ II [Tvi,d^^ , • • •, [Tvi,d^^ ,b] ...] {fug) I 


with constants depending upon how small the aperture of the chosen cones is. 
Notice that the test function / := fufr has the Fourier support as required 
in part (2) oi the statement of Lemma 4. 


Now build cones Dg from the and the other chosen cones as well 
as operators TgUs- Notice that the commutators • • •, [Tvl,Dy^,b] ...] 

and [Ti Di, • • •, [Ti,Di, ^] • • •] reduce signihcantly when applied to a test func¬ 
tion / with Fourier support like ours. When the operators or any 

tensor product Tg Os directly on /, the contribution is zero due to op¬ 
posing Fourier supports of the test function and the symbols of the opera¬ 
tors. The only terms left in the commutators ..., • •](/) and 

, • • •, [Tvl,D.^ ,b]...] if) have the form 0^ (bf) and (X), (bf) 

respectively. 


By repeated use of Lemma 5 we have the operator norm estimates for any 
symbol 6 , valid on the subspace of functions with Fourier support as described 
for /: II Tj^_D^ 6 || 2^2 = II (X)^ T^s,Dvsm2-^2- We conclude that a normalised test 
function / with Fourier support as described in the statement (2) of Lemma 4 
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exists, so that || Ts^Ds{bf )\\2 ^ ||&||BMOi(K'i)- particular, we get the desired 
estimate in (1). QED 

It does not seem possible to pass directly to a lower commutator estimate for 
tensor products of Riesz transforms from that for tensor products of cone op¬ 
erators. Just using tensor products of operators adapted to cones merely gives 
us some lower bound where we are unable to control that a Riesz transform 
does appear in every variable such as required in (3) oi Theorem 4. The reason 
for this will become clear as we advance in the argument. Instead of using op¬ 
erators Tg directly, we will build upon them more general multi-parameter 
Journe type cone operators not of tensor product type that we now describe. 

Let us explain the multiplier we need for i copies of when all dimensions 
are the same. We will explain how to pass to the case of i copies of varying 
dimension dk below. A picture illustrating a base case, a product of two 1- 
spheres, can be found at the end of this section. 

For 0 < b < a < 1, let ip : [—1,1] ^ [—1,1] be a smooth function with 
ip{x) = 1 when a < x < 1 and ip{x) = 0 when b ^ x ^ 0. And let ip be 
odd, meaning antisymmetric about t = 0. The function p> gives rise to a zonal 
function with pole on the hrst copy of denoted by This is 

the multiplier of a one-parameter Calderon-Zygmund operator adapted to a 
cone U(^i, r) for r = 7r/2(l — a). For i > 1 we dehne Ck{Ci^ ■ ■ ■ ■ ■ ■ ■> Vk) 

ioT 1 < k ^ i inductively. In what follows, expectation is taken with respect 
to traces of surface measure. When r]i = +^j, then conditional expectation is 
over a one-point set. 

Cki^i ,... • • •)%) 

(C'fc—1 ) • • • 5 1) hi) • • • ) dk—l) I , ^k—l) d(p]kj Cfc)) • 

If the dimensions are not equal take d = and imbed into 

by the map ^ ■ ■ ■ y^dj) ('^ij • • • > 0,..., 0). Obtain in this manner 

the function Ci and then restrict to the original number of variables when the 
dimension is smaller than d. 

The multiplier J = C'i(^; ■) gives rise to a multi-parameter Calderon-Zygmund 
operator of convolution type (but not of tensor product type), Tj = 

In fact, it is dehned through principal value convolution against a kernel Kj = 
Kci{^..){xi,.. .,Xi) such that 

V/ : I Kj{xi ,..., Xi)dxi = 0, VO < a < /3, Xj e hxed Vj A I, 

Ja<\xi\<0 




1 


1 



-Kj{xi, ... ,Xi)| < A,^|xl| 


\x, 


I di 71 i 


rij ^ 0 . 
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This kind of operator is a special case of the more general, non-convolntion 
type discussed in Section 5. It has many other nice features that will facilitate 
our passage to Riesz transforms. One of them is its very special representa¬ 
tion in terms of homogeneous polynomials, the other one a lower commutator 
estimate in terms of the BMOj norm. 

Lemma 6 Let Ci he a multiplier in as described above, with any 

fixed direction and aperture. Let m he an integer of order d = maxdk- For 
any 5 > 0, the function Ci has an approximation by a polynomial Cfi in the 
nLi4 variables {Ylkn^k^iVkjk | 1 ^ jfc < 4} so that 110* - Of ||cm(§dfc-i) < S 
in each variable separately. 

Qm in^Jexes the norm of uniform convergence on functions that are m times 
continuously differentiable. On the space side. Of corresponds to an operator 
that is a polynomial in Riesz transforms of the variables 

Lemma 7 We are in with partition X = Let T consist of vectors 

^ = (4)Li with fk e Tfc. Let consist of = (4)fcs/,- Let us consider 
the class of Journe type cone multipliers Js = •) of aperture with 

associated multi-parameter Calderon-Zygmund operators Ts^j„. Then we have 
the two-sided estimate 

WHbMOxO^'^) ~ sup ||[Ti^j^, . . . , 6] . . .]|f 2 (Rd)Q < ||&||nMOi(Rrf)- 


In order to proceed with the proof of these lemmas, we will use some well 
known facts about zonal harmonics. Fix a pole f e The zonal harmonic 
with pole ^ of degree n is written as Z^\ri). With t = (fi^rj) e [—1,1], one 
writes Z^'^\r]) = Pn{t) where Pn is the Legendre polynomial of degree n. It is 
common to suppress the dependence on d in the notation for and Pn- 

Z^'^^ are reproducing for spherical harmonics of degree n, _ When 
is harmonic and homogeneous of degree n with = 1 and Y^^\Rr]) = 

YC)[rfi for any rotation R e 0{d) with R^ = then = Z^\ 

The lemma below will aid us in understanding the special form of the functions 

C,. 

Lemma 8 Let fi,f 2 ^ We have 

WMWM = I = 4(2,V2)) 

-E„(Z£>(a2) I 46,«*) = £(({,.1,0)- 


18 


Proof. The first equality is a change of variable, thanks to symmetry of the 
zonal harmonic in its variables and invariance with respect to action of the 
measure preserving elements of the orthogonal group hxing poles or that 
we now detail. By a rotation in one of the spheres, assume = ^2 = f - Take 
a small ball 

^$,,, 1 ( 02 ; ^ 2 ) = {02 : d{a 2 ,al) < 62 } n {oa : d{a 2,0 = ^( 771 ,01- 

Note {oa : d(a2,0 = Every Oa e gives rise to a 

canonical orthogonal map cxaj along geodesics in a scaled copy of Lifted 
to these are orthogonal maps hxing 0 Let cr° hx ^ and map 0 ° to rji. Let 
a? = o-°(h 2 )- We observe that {a°aa 2 (h 2 ) : 02 e £ 2 )} = S^,^ 2 (a?;£i) 

with El so that 

P(d(a2,a^) < ea I c ?(0 «2) = c?(Ohi)) = P(d(ai,a?) < ei | d{f,ai) = d( 0 h 2 ))- 

Together with the symmetry and the rotation property Z^\a) = Z^"‘\r]) = 
Z^l^{a{r])), we obtain the hrst equality. 

For hxed Oi, the function Z^f\ai) = ^(”^( 771 ) is a function harmonic in Mf, 
TT-homogeneous. These properties are preserved when taking expectation in 
oi. So the expression E(Z^”)(ai) | d(^i,ai) = ^(^ 2 ,^ 2 )) remains harmonic 
(regarded as a function in R"^), 7 r-homogeneous. From the form E,{Zlj^^{a 2 ) \ 
^( 1 ^ 2 , 02 ) = d(.^i, 77 i)) we learn that its restriction to depends only upon 
d(^i, 771 ). This implies that it is a constant multiple of the zonal harmonic with 
pole ^ 1 . Exchanging the roles of 771 and 772 gives 

E(Z("l(ai) I dKi.Oi) - d(( 2 ,ri 2 )) - c„d”Ti)Zi 7 '(>; 2 )- 

When assuming the normalization — 1 then c„ — 1. 

This is a gernalisation of the classical symmetrising of the cosinus sum formula 
l/2(cos(a; + y) + cos(x — y)) = cos(a;) cos{y). 


QED 

Proof, (of Lemma 6 ) It is well known that zonal harmonic series have conver¬ 
gence properties when representing smooth zonal functions similar to that of 
the Fourier transform. For any given m and sufficiently smooth (p of the type 
described above, then 

<^ 1 ( 6 ; hi) 

n 

where the convergence is C^-uniform. The degree of smoothness required for 
ip to obtain convergence in the C"* in the above expression depends upon m 
and the dimension d. For our purpose, we choose m ^ d. 
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Let us denote this function’s representation of degree by a series of zonal 
harmonics by ??i)- 

n^N 

For every 5 > 0 there exists N so that we have the estimate 

l|cr>(6;>)i)-C’i(6;>?i)llc»(s'.->) <^- 

In the case of i copies of spheres, we define inductively in the same 

manner as C*. Let us for the moment make all dimensions equal using the 
argument discussed above. So we set 

= Eafc_i(C'S(^i,.. .,ak-i]Vi, ■ ■ -^Vk-i) I d(afc_i,^fc_i) = d{r]k,^k))- 

We claim the identity 

'dnflzl'^\r]k). ( 2 ) 

n^N k^l 


This is trivially true for i = 1. For i > 1 induct on the number of parameters: 

= Eai_i ( Yj ^nY\^c^\vk) I d{ai-i,^i-i) = d{r]i,^i) 

\n^N k=l 
i-2 

n^N k^l 

= Y 

n^N k=l 


The hrst equality is the dehnition of the second one is the induction 

hypothesis and the last an application of Lemma 8. 

It follows that neither Ci nor depend on the order chosen in their dehni¬ 
tion and 

n fc=l 

where the convergence is in in each variable. 
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Next, we study the terms arising in multipliers of the form When all 

dimensions are equal, indeed, nh. 4 ”’ {rjk) has the important property that, 
as a product of n-homogeneous polynomials, has only terms of the form 

k^l k^l \ife = l / 

where rjk e and ak = {c^kj,,) are multi-indices with la^l = = n 

for all k. This form is inherited by with varying n. It shows that is 
indeed a polynomial in the variables OI=i Vk,jk- When the dimensions dk are 
not equal, observe that by restricting back to the original number of variables, 
we certainly lose harmonicity of the polynomials, but not n-homogeneity or 
the required form of our polynomials. QED 

Proof, (of Lemma 7) By Lemma 4 we know that for each parameter 1 4 
s 4 I there exists a tensor product of cones Dg = (S)keis DiCk,rk) with : = 
Sfcs/s ^ £ Tf'fc and test functions fs supported as described in 

Lemma 4 part (2) so that 

||[rL,Di, . . . , [Ti^Di,b] . . .](/)||2 ^ ||^||BMOx(R‘i) 

where / = fs- We make a remark about the apertures r^. Let 

denote geodesic distance on where antipodal points have distance tt. Let 
be the set of directions of Dg. Remember that according to Lemma 4, 
one component had a specihc direction e and possibly large aperture 
with (1 -I- < 7 r/ 2 . Let us choose the other directions arbitrarily but with 

apertures small enough so that (l-fr)(r(,^^ -I-(z —< n/2. Now choose 

an aperture Vg < t^I2 so that (1 -I- < 4 < 7 r/ 2 . 

Writing ig = l/^l, we hnd Journe type cone multipliers Jg = ac¬ 

cording to the construction above with center and aperture Vg. 

We are going to observe that = 1 on spt(Ds) and = — 1 on the Fourier 
support of fg. Let us drop the dependence on s for the moment. We see in an 
inductive manner that Cj(^; ■) takes the value 1 in a certain ball of radius 
r < ti/2 centered at We show that 

Y^d{fk,r]k) < r ^ C'i(^, 77 i,..., 77 i) = 1 . 

k 

When z = 1, the assertion is obviously true: d{fi, rji) < r ^ hi) = 1 by 

the choice of (p,r and dehnition of Ci. For z > 1, we proceed by induction. 
Assume that I]I=\< r implies ..., 771 ,..., = 1 . 

Let us assume that Yjk=i d{fk, Vk) < r. Remembering the dehnition of C'j(^; •) 
we assume d{ai-i,fi-i) = d{rii,^i). By the triangle inequality Xil=ihfc) + 
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d{ai-i,r]i_i) < 'Zl=i(^i^k,Vk)+d{ai^u^i-i)+d{^i_i,r]i_i) = Yjk^i^i^k^Vk) < r. 
So 

Ci-i(^i, ^2; • • •! Oj_i; r^i,..., rji^i) = 1 

for all aj_i relevant to the conditional expectation in the dehnition of C'j(^; ■). 
The statement for i follows. 

Since C'j(^; ■) does not depend npon the order of the variables in its constrnc- 
tion, we are also able to see exactly as done above that when ak = —I for 
exactly one choice of k, then d{akiki Vk) < r ^ C'i(^; 771 ,... , r/,) = -1. 

Consider associated mnlti-parameter Calderon-Zygmnnd operators and 
Ids = (S)keis identity on the variable. Now 


..., [Ti,j„b] ...](/) = + Id ,,..., [Ti^j, + Idi, b] ...](/) 

= ®{T,J. + Id,)(h!) 

s=l 

With II (S)s^,{Ts,j^ + Ids){bf)\\2 ^ II (S)s^,Ts,nAM)\\2 and 0 [^,Ts,dM) = 
[^1,01) • • •) [Ti,Di, b] ...](/) we get the desired lower bonnd on the Jonrne com- 
mntator as claimed. QED 


Let ns illnstrate the base case (S^)^ by a pictnre. The pictnre is simplihed in 
the sense that the odd fnnction (p above is replaced by an indicator fnnction 
of an interval. 





Cone fnnctions based on the obliqne strips con¬ 
taining ^ are averaged. In the two-dimensional 
case, expectation is over a one or two point 
set only. The rectangle aronnd ^ with sides par¬ 
allel to the axes representing illnstrates the 
snpport of the tensor prodnct of cone opera¬ 
tors with direction The longer side is the 
^ apertnre that arises from the Hankel part. The 
short sides can be chosen freely as they arise 
from the Toeplitz part and are chosen small so 
that the rectangle fits into the obliqne sqnare. 
The other small rectangle corresponds to the 
Fonrier snpport of the test fnnction /. 


Proof, (of Theorem 4) 

In contrast to the Hilbert transform case, both lower bonnds reqnire separate 
proofs. This is a notable difference that stems from the loss of orthogonal 
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subspaces in conjunction with the special form of the Hilbert transform only 
seen in one variable. It does not seem possible to get a lower estimate (3)^(2) 
directly. 

(1) ^(2). The upper bound (1)^(2) is an easy consequence of the upper 
estimates of iterated commutators of single Riesz transforms. The lower bound 

(2) ^(1) follows from a standard fact on multipliers in combination with the 
main result in [18], the two-sided estimate for iterated commutators with Riesz 
transforms, similar to the hrst arguments used in 4. 

(1 )<^ (3). The upper bound (1 (3) follows from the tensor product structure 
and use of the little product BMO norm (see also the remarks in section 7). 
The lower bound (3)^(1) uses the considerations leading up to this proof: 
Suppressing again the dependence on s, we see that the multiplier Q is an 
odd, smooth, bounded function in each rjk when the other variables are fixed. 
Furthermore, since (p, written as a function of t rj) is odd with respect 
to f = 0, the above series has ipn ^ 0 at most when n is odd and thus is 
odd. So is as a sum of odd functions odd. 

We are now also ready to see that Tj, the Journe operator associated to 
the cone J = ■) as well as the operator associated to are 

paraproduct free. In fact, applied to a test function / = (x)^ with acting 
on the variable and where /z = 1 for some I gives Tj{f) = 0. To see this, 
apply the multiplier ■) in the I variable (acting on 1) hrst, leaving the 

other Fourier variables hxed. The multiplier function 

r]i ^ (I; hi, • • •, ^i) = 2 fl 

n^N k^l,k^l 

is, as a sum of odd functions, odd on bounded by 1 and uniformly 

smooth for all choices of rjk with k 1. As such it gives rise to a paraproduct 
free convolution type Calderon-Zygmund operator in the Zth variable whose 
values are multi-parameter multiplier operators. 

Due to the convergence properties proved above, the difference 

>(?;.) 

gives rise to a paraproduct free Journe operator with Calderon-Zygmund norm 
depending on N. This is seen by an application of an appropriate version of 
the Marcinkievicz multiplier theorem. 

By our stability result on Journe commutators in section 5, Corollary 3, there 
exist for all 1 ^ s < / integers Ng so that with e is a 

characterizing set of operators via commutators for BMOi(M'^). This is a hnite 
set of possibilities because of the universal choice of the r* and hniteness of 
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the set T. Using the multi-parameter analog of the observation \AB^h\ = 
A{B^h\ + {A^h\B and the special form of the ■), leaves us with the 

desired lower bound: Observe that when [AS, h\ has large norm then either 
[A, h\ or [S, h\ has a fair share of the norm. We use this argument hnitely many 
times in a row for operators that are polynomials in tensor products of Riesz 
transforms This hnishes (3)^(1). QED 

We remark that there are two cases of dimension greater than 1, where the 
proof simplihes. In the case of arbitrarily many copies of one can work with 
the multiplicative structure of complex numbers and avoid the symmetrizing 
procedure to obtain cone functions with the appropriate polynomial approxi¬ 
mations. If the dimensions are arbitrary, but only tensor products of two Riesz 
transforms arise, one can avoid part of the construction above by using the 
addition formula for zonal harmonics. 


5 Real variables: upper bounds 


In this section, we are interested in upper bounds for commutator norms by 
means of little product BMO norms of the symbol. In the case of the Hilbert 
transform, we have seen that these estimates, even in the iterated case, are 
straightforward. Other streamlined proofs exist for Hilbert or Riesz trans¬ 
forms when considering dyadic shifts of complexity one, see [29], [30] and [19]. 
When considering more general Calderon-Zygmund operators, the arguments 
required are more difficult, in each case. The hrst classical upper bound goes 
back to [6], where an estimate for one-parameter commutators with convolu¬ 
tion type Calderdn-Zygmund operators is given. Next, the text [18] includes 
a technical estimate for the multi-parameter case for such Calderdn-Zygmund 
operators with a high enough degree of smoothness. This smoothness assump¬ 
tion was removed in [8] thanks to an approach using the representation formula 
for Calderdn-Zygmund operators by means of inhnite complexity dyadic shifts 
[15]. This last proof also gives a control on the norm of the commutators which 
depends on the Calderdn-Zygmund norm of the operators themselves, a fact 
we will employ later. Below, we give an estimate by little product BMO when 
the Calderdn-Zygmund operators are of Journe type and cannot be written as 
a tensor product. While this estimate is interesting in its own right, remem¬ 
ber that it is also essential for our characterization result, the lower estimate, 
in section 4. The first generation of multi-parameter singular integrals that 
are not of tensor product type goes back to Fefferman [11] and was general¬ 
ized by Journe in [16] to the non-convolution type in the framework of his 
T(l) theorem in this setting. Much later, Journe’s T(l) theorem was revis¬ 
ited, for example in [22], [27], [28]. See also [23] for some difficulties related 
to this subject. The references [22] in the bi-parameter case and [28] in the 
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general multi-parameter case include a representation formula by means of 
adapted, infinite complexity dyadic shifts. While these representation formu¬ 
lae look complicated, they have a feature very useful to us. ‘Locally’, in a 
dyadic sense, they look as if they were of tensor product type, a feature we 
will exploit in the argument below. We start with the simplest bi-parameter 
case with no iterations and make comments about the generalization at the 
end of this section. 

The class of bi-parameter singular integral operators treated in this section is 
that of any paraproduct free Journe type operator (not necessarily a tensor 
product and not necessarily of convolution type) satisfying a certain weak 
boundedness property, which we dehne as follows: 

Definition 2 A continuous linear mapping T : C®(M”')(x)C'®(R™') [C'®(R"')(x) 

C®(R"*)]' is called a paraproduct free bi-parameter Calderon-Zygmund oper¬ 
ator if the following conditions are satisfied: 

1. T is a Journe type bi-parameter 5-singular integral operator, i.e. there exists 
a pair {Ki, K 2 ) of 5CZ-5-standard kernels so that, for all fi, gi e (^“(R”) and 

f2,g2^C^{^n, 

(T{fi^ f2),gi® 92 ) = J fi{yi)(Ki{xi,yi)f2,g2}gi{xi)dxidyi 

when sptfi n sptgi = 0; 

(T{fi0 f2),gi0 92 ) = J f2{y2KK2{x2,y2)fi,gi)92{x2)dx2dy2 

when sptf 2 n sptg 2 = 0. 

2. T satisfies the weak boundedness property \(T{xi 0 Xj)^Xi ® Xj}\ ~ 
for any cubes J c R"', J e R™. 

3. T is paraproduct free in the sense that r(l (x) •) = T(- (x) 1) = T*(l (x) ■) = 
T*(-(x)l) = 0. 


Recall that a hCZ-^-standard kernel is a vector valued standard kernel taking 
values in the Banach space consisting of all Calderon-Zygmund operators. It 
is easy to see that an operator dehned as above satishes all the characterizing 
conditions in Martikainen [22], hence is bounded and can be represented 
as an average of bi-parameter dyadic shift operators together with dyadic 
paraproducts. Moreover, since T is paraproduct free, one can conclude from 
observing the proof of Martikainen’s theorem, that all the dyadic shifts in the 
representation are cancellative. 

The base case from which we pass to the general case below, is the following: 
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Theorem 6 Let T be a paraproduct free bi-parameter Calder on-Zygmund op¬ 
erator, and b be a little bmo function, there holds 

II I|l2(R"xR'") 0 ~ ll^ll 6mo(K'»xK"*)) 

where the underlying constant depends only on the characterizing constants of 

T. 

Proof. 

According to the discussion above, for any sufficiently nice functions f, g, one 
has the following representation: 

00 00 

(Tf,g} = 2 2 (3) 

1202=0 

where expectation is with respect to a certain parameter of the dyadic grids. 
The dyadic shifts are defined as 

^nii*2j2 j 

- S S S S 0‘hJiKxi2J2K2{f 1 hh ® hif}hj^ (g) hj^ 

Ki£T>i K2e'D2 I2,J2'^K2,I2,J2^T^2 

l(h)=2-nt{Ki) i(l 2 )^ 2 -i 2 i{K 2 ) 

e{Ji)=2-ne{Ki) e{J2)=2-n£(K2) 

(hdl) (*2j2) 

-suss (^iiJiKii2J2K2(f^ hh ® hj^yhj^ (g) hjj. 

Kl h,Jic^Ki K2 l2,J2^K2 


The coefficients above satisfy ajj^j.^Kii2J2K2 ^ which also guaran¬ 
tees the normalization ^ 1 - Moreover, since T is paraproduct 

free, all the Haar functions appearing above are cancellative. 

It thus suffices to show that for any dyadic grids Pi, V2 and hxed ii,ji, *2,^2 ^ 
N, one has 

||[6, < (1-h max(A,ii))(l + max(z 2 , j 2 ))||&||bmo||/||L 2 , (4) 

as the decay factor Ji)^ 2-™a,xp2,i2) ('3^ guarantee the convergence 

of the series. 

To see ( 4 ), one decomposes h and a Lf test function / using Haar bases: 

[&, =2] <6, hi, (g) hi,Xf, hj, (g) hj.'yihi, (g) hi„ (g) hj,. 

h,l2 Jl,J2 


A similar argument to that in [8] implies that [hj, (g) hi^, (g) hjj is 
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nonzero only if Ii c or I2 c where denotes the Zi-th dyadic an¬ 
cestor of Ji, similarly for Hence, the snm can be decomposed into three 
parts: Ji c and I2 c (regnlar), Ji c and I2 2 J2‘^\ h 2 
and I2 c (mixed). 


1 ) Regular case: 

Following [8] one can decompose the arising snm into sums of classical bi¬ 
parameter dyadic paraproducts B^ip, f) and its slightly revised version /): 

for any integers A:, / ^ 0, B^ i is the bi-parameter dyadic paraproduct dehned 
as 


BkAbJ) = 2/3/X^ V) ® ® ® 

I,J 


where ( 3 jj is a sequence satisfying \f 3 ij\ ^ 1 . When k > 0 , all Haar functions 
in the hrst variable are cancellative, while when k = 0 , there is at most one 
of being noncancellative. The same assumption goes for the second 

variable. Observe that when k = I = 0 , B^ i becomes the classical paraproduct 
Bq. It is proved in [8] that 




with a constant independent of k, I and the product BMO norm on the right 
hand side. 


Then since little bmo functions are contained in product BMO, this part can 
be controlled. More specihcally, write 


[b, = 2 Z <^’ ® hi,Xf, hj, ® hjXhr, (X) ® hj,) 

h,l2 Jl,J2 

- S S ® hiXif, hj, (X) ® hi.hj,) 

h,l2 Jl,J2 

= : / + //, 


then one can estimate term I and II separately. According to the dehnition of 
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dyadic shifts, term I is equal to 


(b, hi^ (X) hj^Xf, ^Ji ® hjXhi^ (X) hj^- 


( S S 

£(J()=2n-Jl£(Ji) £(J')=2»2-J2£(J2) 




hri (S)h 


(il) ( 12 ) 

- s s s s s (b, hi^ (X) hiXif, hj^ (X) hj 2 >h/i (x) 

Jl-.Jl<^Ki J 2 -.J 2 ^K 2 IvhcKi I 2 .I 2 ^K 2 


(il) (i2) 

( 2j 2j (^JlJ[K^J2J^K2hj[®hj!^j 

J[-.J[czKi J':J'cX 2 

(*i,ii) ( 12 , 72 ) 

= Y,{bXh®hiXhi,®hj,Y, S S «JiJ(XiJ 2J'K2(/5 ^Ji ® hjXhj!^ (X) /i 

q,/2 Ji,J'cKi J2 ,J'cK2 

i^2^/2 

= 2 ® ^i2)hh ®hi^ Yj Zj ® ® ^J'• 


Because of the supports of Haar functions, the inner sum above can be further 
decomposed into four parts, where 


Il,l2 J'2^h h,l2 J[^Il J': J'c/2C 

‘"-n s n- s i; 


Hence, using the same technique as in [8], one has 

/ = 2 hi, (8) h\ (8) h\)hi, (8) h,2|/lr'/'|/2r'/^ 

hj 2 


which is a bi-parameter paraproduct BqXX)- Moreover, one has 
II = Y(b,hi,®hjXhj,^hi, Y <S^^^^'^'XX],®hjX\Ii\-h^hj, 

jpj'c 72 C=j'^^ 2 ) 

= Z Z ® h], 8) (8) hXh\-h^\jY\-h^ 

1=0 h A 

= YBoi{b, ^™V), 

1=0 

where constants / 3 j' e {1,-1}, and Bqi are the generalized bi-parameter para- 
products of type (0, /) dehned in [8] whose operator norm is uniformly 


bounded by ||6||bmo product BMO. Similarly, one can show that 
III = 2 

k =0 k= 0 l =0 


Since II^IIbmo ^ l^lbmo, ah the forms above are bounded. This completes 
the discussion of term I. 

To get an estimate of term II, we need to decompose it into hnite linear combi¬ 
nations of /)). By linearity, one can write on the outside 

from the beginning, and we will only look at the inside sum. One splits for 
example the sum regarding the hrst variable into three parts: Ii ^ Ji, Ii = Ji, 
Ji c c If we split the second variable as well, there are nine mixed 
parts, and it’s not hard to show that each of them can be represented as a 
hnite sum of Bki{b, /). We omit the details. 


2 ) Mixed case. Let’s call the second and the third ‘mixed’ parts, and as the two 
are symmetric, it suffices to look at the second one, i.e. Ii c 
In the hrst variable, we still have the old case Ji c that appeared in 
[8], so morally speaking, we only need to nicely play around with the stronger 
little bmo norm to handle the second variable. For any hxed /i, Ji, I2, J2, since 
I2 2 J2^\ the dehnition of dyadic shifts implies that 

hi, (X) (X) hjJ = (X) hi,hj,) 


and 

(X) hi,hj,) = (X) hjJ. 

Hence, we still have cancellation in the second variable, which converts the 
mixed case to 

T T (b, hi, ® hi,Xf, hj, (X) hjX[hi„S^^^^^^^^]{hj, (X) hM 

hczjX^ 

- T hj, » hjIAh,,, (hj, 0 <6. h,, ® h) (.i2))h^j(i2)h.h) 

= hi, ® hj, ® hjX[hi„S^^^^^^^-]{hj, ® hj,) 

= S Y/^{b)jip)MXi{f. hj, (X) hjX{hi„S^^^^^^^%hj, (X) hjJ, 
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where (h) M2) denotes the average value of b on J2^\ which is a function of 
'^2 

only the hrst variable. 


In the following, we will once again estimate the hrst term and second term of 
the commutator separately, and the norm of each of them will be proved 
to be bounded by ||6||bmo||/||L2. 

a) First term. 


By dehnition of the dyadic shift, the hrst term is equal to 


j(i2); ^Ji ® hj^}- 

e{j'^)=2n-h£{j^) e{j^)=2^2-32e(j2) 




which by reindexing Ki := is the same as 


Il,J 2 

(*i) (ii) O2) 

Ky.Ki^h 

= 2«&V2),/iq>i/iq 2] hj,®{S^^^^^^^^{{fAj2')2®hj,)Aj[')u 


h,J 2 




where the inner sum is the orthogonal projection of the image of (/, hj^')2®hj^ 
under onto the span of {hj/} such that 3 Ii. Taking into account 

the supports of the Haar functions in the hrst variable, one can further split 
the sum into two parts where 

"-- 1 . S ■ 

J2 Il^J[ J2 

Summing over Jj hrst implies that 

‘ '*'.>1'*'. ® ® <S""‘””«/. hjM- ® hj,),h]M) 

J2 h ' 

= : Bo«&>^(«2), hj,}2 ® hj,)) 


where B^ib^ f) := 2/(^5 a classical one-parameter para- 

product in the hrst variable. Note that its norm is bounded by ||b||BMo||/||L2. 
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Moreover, according to the definition of for any fixed J2 

Snnnn^(^f^ hj,>2 (x) hj,) = ^ (g) 

If. r'C'2)_ t(»2) 

'■'2 

In other words, hj^^^hj^) only lives on the span of {hj/ : = 

J2*^^}. Hence, by linearity there holds 

I = Y, 2 i?0(<fe>^(«2)X^™^'H</,/ij2>2® 

J2 jf. T'(h)_ Ah) 

'^2"^2 '^2 

■^2 J2:4*2)^j'(J2) 

Thus, orthogonality in the second variable implies that 

II tll^ 

II-' llL2(RnxRm) 

4 J2:4*2)^j'(J2) 

~ y~! iK^) r'(J2) IIbMO(R’*') ^ (/) ^72)2 ® hj2)) ^J^)2||i2(Rn). 

4 J2:4*2)^j'02) 

Observing that |Kb>^/02) ||bmo(R") < <||fe||BMO(R")>j'02) < l^lbmo, one has 

< IlftllLoIllKS’”™^ H </.A*>2®AA).ftj;>2lli.(.», 

4 J2;4*2)^j'02) 

= ll^llbmoll ^ (/; ^J 2)2 ® hj 2 )) ^J^)2 ® hj' ||^ 2 (RnxRm). 

■^2 J2:4*2)^j'02) 

Note that the sum in the norm is in fact very simple: 

4 J2:4*2)^j'02) 

-H H <S‘‘"“«/.fc*>2®A*),A2;>2 0Aji 

J2 7/. t'0’2)_ t(^ 2) 

’^2''^2 ~^2 

= 2 hj,>2 ® hjJ = 

J2 

Hence, the uniform boundedness of the operator norm of dyadic 

shifts implies that 

I|-^IIl2(R’*xR’") ~ ll^llbmoll/llL2(RnxR™)- 

In order to handle H, we split it into a finite sum depending on the levels of 
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1 1 upon J(, which leads to 


II = 2 /ij 2>2 (X) hjJ, hj.>i 


fc=0 J 2 J{ 
ji 


hj,}2 ® hjJ, hjj>l 

fc=0 J2 J( 


fc=0 J 2 


where Bk{b, f) := J]i I 3 i,k(b, hj(k)Xf, hi}hi\X'‘^\ is a generalized one-parameter 
paraproduct studied in [8], whose norm is uniformly bounded by ||&||bmo||/||l 2 , 
independent of k and the coefficients j 3 i^k ^ {I; ~ 1 }- Then one can proceed as 
in part I to conclude that 

||-^-^||L2(RnxRm) < ( 1 -t" jl ) || 6 1 | bmo || / || L2 (K’^ X K"* ) ) 
which together with the estimate for part I implies that 


First term||i2('Rn><Rm) < (1 + Ji)||fe||bmo||/||L2(R"x]R"')- 


b) Second term. 

As the second term by linearity is the same as 

Qiljli2j2 (S S ((b)j(i 2 ),hh)i(f, hji ® hjXhiXji ® , 

■^2 

the boundedness of the shift implies that it suffices to estimate the 

norm of the term inside the parentheses. Since /i n A ^ 0, one can further 
split the sum into two parts: 

J2 /l£Jl J2 

Summing over Ji first implies that 

I = J]Y2^(I}j(i2),hiXi(f, h]^ 0 hjXhi.h], 0 hj^ 

J2 h ' 

= : 2]^o«&> 02),</, hj2>2) 0hj2, 

J 2 ^ 

where Bo{b, f) ■= is a classical one-parameter para- 
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product in the first variable. Hence, 


II tll^ 

F llL2(]Rnx]Rm) 




J2 


< 


\\(.t)}j(i2) ||bMO(]R")IK/) ^J2)2|Il2(R" 


J2 


< 


II ^llbmo SiK/ ! ^J2)2||i2(Rr>) = ll^ll bmo ll/lli2(R’i 


xR"*)- 


J 2 


For part II, note that it can be decomposed as 


h 

II = j(»2)5 hjw')i'\fi ® hj^h (X) hj^ 

/c —0 J2 Ji 

hjA2, hjAihj, (X) hj, 

/c —0 J2 Ji 

= ■ 2 ] 2 ]^fc«&>j(i 2 ),</>^J 2 > 2 )®^J 2 > 

Ai—0 J2 

where coefficients ^ {I5 ~1} and the norm of the generalized para- 
product Bk is uniformly bounded as mentioned before. Therefore, the same 
argument as for part I shows that 


11 -f.f ||l2(R"xR’") ^ (1 + H)||&||bmo||/||L2(]RnxR™), 

which completes the discussion of the second term, and thus proves that the 
mixed case is bounded. QED 


The upper bound result we just proved can be extended to to arbitrarily 
many parameters and an arbitrary number of iterates in the commutator. 
To do this, consider multi-parameter singular integral operators studied in 
[ 28 ], which satisfy a weak boundedness property and are paraproduct free, 
meaning that any partial adjoint of T is zero if acting on some tensor product 
of functions with one of the components being 1. And consider a little product 
BMO function b e BMOi(M‘^). One can then prove 

Theorem 7 Let us consider d = (di, ... ,dt) with a partition X = (A)i<jsisz 
of { 1 ,... ,t} as discussed before. Let b e BMOx(RA '^s denote a multi¬ 

parameter paraproduct free Journe operator acting on functions defined on 
( 8 )fce/s have the estimate below 

m,...[Ti,b]...]\\ L2(Rd)Q < ||fe||BMOi(Rrf)- 


The part of the proof that targets the Journe operators proceeds exactly the 
same as the bi-parameter case with the multi-parameter version of the rep- 
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resentation theorem proven in [ 28 ]. Certainly, as the number of parameters 
increases, more mixed cases will appear. However, if one follows the corre¬ 
sponding argument above for each variable in each case, it is not hard to 
check that eventually, the boundedness of the arising paraproducts is implied 
exactly by the little product BMO norm of the symbol. The difficulty of higher 
iterates is overcome in observing that the commutator splits into commutators 
with no iterates, as was done in [8]. We omit the details. 

The assumption that the operators be paraproduct free is sufficient for our 
lower estimate. The general case is currently under investigation by one of the 
authors. Important to our arguments for lower bounds with Riesz transforms 
is the corollary below, which follows from the control on the norm of the 
estimate in Theorem 7 by an application of triangle inequality. It is a stability 
result for characterizing families of Journe operators. 

Corollary 3 Let for every 1 ^ s ^ I be given a collection Ts = {Tgj^} of 
paraproduct free Journe operators on that characterize 

via a two-sided commutator estimate 

\\H BMOxhS.’i) ~ ■ ■ •]llL2(R<i)Q ^ \\b\\ BMOx{R<i) ■ 

3 

Then there exists e > 0 such that for any family of paraproduct free Journe 
operators T's = {TfjJ with characterizing constants lir^jJIcz ^ £? the family 
TfjJ still characterizes 


6 Weak Factorization 


It is well known, that theorems of this form have an equivalent formulation 
in the language of weak factorization of Hardy spaces. We treat the model 
case and BMO(i 3)(2)(M‘^) only for sake of simplicity. The other 

statements are an obvious generalization. For the corresponding collections of 
Riesz transforms Tlk,jk ^ ^ BMO(i3)(2) (M^), 1 < s < 3 , by unwinding the 
commutator one can define the operator Hj such that 

<[R2 ,j2 , b]]f, 9 )l^ = (b, ^jif, 9 )}l^- 


Consider the Banach space of all functions in of the form 

/ = normed by 

||/||l2*l 2 = in ({SHidibidiw 

3 * 

with the inhmum running over all possible decompositions of /. Applying a 
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duality argument and the two-sided estimate in Corollary 2 we are going to 
prove the following weak factorization theorem. 

Theorem 8 ^ (g) coincides with 

the space . In other words, for any f e + 

there exist sequences (fiji/jf e such that f = 2]^ X!i V'/) 

and ll/ll - SjSi Il</>il|2||^f||2. 

Proof. Let’s hrst show that is a subspace of (g) + 

L^(M'’*i)(x)if;^g(MC2,rf3)^_ Recalling the remark after Theorem 2 , this is the same 
as to show V/ e / is a bounded linear functional on BMO(i3)(2) (M^). 

This follows from the upper bound on the commutators since 

(b, g g; TIM, ii)) - [RbMTh. iM.4}. 

j i j i 

Now we are going to show 

f {l J/^h \\f\\L^*L^ ^ l} II^I|bMO(i 3 )( 2 ) 

which gives the equivalence of (g -f L^(M'^i)(x)iL^g(MC2.c«3)^ 

norm and the norm, thus showing that the two spaces are the same. 

To see this, note that the direction < is trivial, and the direction > is implied 
by the lower bound of commutators. For any b e BMO(i3)(2), there exists j 
such that ||&||bmO( 13)(2) ~ I|[^2,i2) ^]]||- Hence, there exist e 

with norm 1 such that 

II^IIbmo( 13 )( 2 ) ^ \([R 2 ,j 2 dRi,jiR 3 ,h,b]](j),ij)\ = \(b,nj{(j),ij))\ ^ lhs, 

which completes the proof. QED 


7 Remarks about our results in 


As mentioned before, the two-sided estimates stated in section 4 and in partic¬ 
ular Theorem 5 hold for all 1 < p < oo. The fact that upper estimates hold for 
the Riesz commutator in in the case where no tensor products are present is 
proved in [ 18 ] as well as [ 19 ]. It stems from the fact that endpoint estimates for 
multi-parameter paraproducts hold for all 1 < p < oo [ 24 ], [ 25 ]. This estimate 
carries over easily to tensor products of Riesz transforms or any other tensor 
products of operators for which we have estimates on the commutator: one 
uses [TiT 2,6] = Ti[r2,6] + [Ti, &]T’2 fo handle arising tensor products, followed 
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by a correct use of the little product BMO norm. The argument is left as an 
exercise. 

The lower estimate or the necessity of the BMO condition can be derived 
from interpolation. In fact, suppose we have uniform boundedness of our com¬ 
mutators with operators running through all choices of Riesz transforms and 
some symbol b in L^. Then by duality, we have boundedness in where 
1/p+l/q = 1. In fact, [T,b]*f = —[T*,b]f = — [T*, 6]/ shows that the bound¬ 
edness of adjoints is inherited. The same reasoning holds for iterated commu¬ 
tators of tensor products. Thus by interpolation, the boundedness holds in 
and the symbol function b necessarily belongs to the required BMO class. 


References 

[1] Arlen Brown and Paul Halmos, Algebraic properties of Toeplitz operators, J. 
Reine Angew. Math. 213 (1964), 89-102. t 

[2] Lennart Carleson, A counterexample for measures bounded on spaces for 
the bidisk, Mittag-Leffler Rep. No. 7, Inst. Mittag-Leffler (1974). |5 

[3] Ronald Coifman, Pierre-Louis Lions, Yves Meyer, and Stephen Semmes, Com¬ 
pensated compactness and Hardy space, J. Math. Pures Appl. 72 (1993), no. 9, 
247 -286. 13 

[4] Sun-Yung Chang and Robert Fefferman, Some recent developments in Fourier 
analysis and -theory on product domains. Bull. Amer. Math. Soc. (N.S.) 12 
(1985), no. 1, 1-43.MR 86g:42038 15 

[5] _, A continuous version of duality of with BMO on the bidisc, Ann. 

of Math. (2) 112 (1980), no. 1, 179-201. 15 

[6] Ronald Coifman, Richard Rochberg, and Guido Weiss, Factorization theorems 
for Hardy spaces in several variables, Ann. of Math. (2) 103 (1976), no. 3, 
611-635. T2, 24 

[7] Misha Collar and Cora Sadosky, The Helson-Szegd theorem in of the bidi- 
mensional torus, Contemp. Math. 107 (1990), 19-37. |5 

[8] Laurent Dalenc and Yumeng Ou, Upper bound for multi-parameter iterated 
commutators, arXiv:1401.5994, 1-18. 12, 24, 27, 28, 29, 32, 34 

[9] Laurent Dalenc and Stefanie Petermichl, A lower bound criterion for iterated 
commutators, J. Fund. Anal. 266 (2014), no. 8, 5300-5320. |13, 14 

[10] Robert Fefferman, Bounded mean oscillation on the polydisk, Ann. of Math. 
(2) 110 (1979), no. 2, 395-406. 15 

[11] _, Singular integrals on product domains. Bull. Amer. Math. Soc. (N.S.) 

4 (1981), no. 2, 195-201. 125 


36 





[12] _, Harmonic analysis on product spaces, Ann. of Math. (2) 126 (1987), 

no. 1, 109-130. T5 

[13] Sarah Ferguson and Michael Lacey, A characterization of product BMO by 
commutators, Acta Math. 189 (2002), no. 2, 143-160. |2 

[14] Sarah Ferguson and Cora Sadosky, Characterizations of bounded mean oscil¬ 
lation on the polydisk in terms of Hankel operators and Carleson measures, J. 
Anal. Math. 81 (2000), 239-267. ]2, 3, 5, 6 

[15] Tuomas Hytonen, The sharp weighted bound for general Calder on-Zygmund 
operators, Ann. of Math. 175 (2012), no. 3, 1473-1506. 124 

[16] Jean-Lin Journe, Calderon-Zygmund operators on product spaces. Rev. Mat. 
Iberoamericana 1 (1985), no. 3, 55-91. fd, 25 

[17] _, A covering lemma for product spaces, Proc. Amer. Math. Soc. 96 

(1986), no. 4, 593-598. j 

[18] Michael Lacey, Stefanie Petermichl, Jill Pipher, and Brett Wick, Multiparam¬ 
eter Riesz commutators, Amer. J. Math. 131 (2009), no. 3, 731-769. 12, 3, 5, 
13, 23, 24, 35 

[19] _, Iterated Riesz commutators: a simple proof of boundedness., Contemp. 

Math. 505 (2010), 171-178. 124, 35 

[20] _, Multi-parameter Div-Curl Lemmas, Bull. London Math. Soc. 10 

(2012), 1123 -1131. t3 

[21] Michael Lacey and Erin Terwilleger, Hankel operators in several complex vari¬ 
ables and product BMO, Houston J. Math. 35 (2009), no. 1, 159-183. | 

[22] Henri Martikainen, Representation of bi-parameter singular integrals by dyadic 
operators, Adv. Math. 229 (2012), no. 3, 1734-1761. 14, 25 

[23] Henri Martikainen and Tuomas Orponen, Some obstacles in characterizing the 
boundedness of bi-parameter singular integrals, arXiv: 1404.2216, 1-10. |25 

[24] Camil Muscalu, Jill Pipher, Terence Tao, and Christoph Thiele, Bi-parameter 
paraproducts, Acta Math. 193 (2004), no. 2, 269-296. 135 

[25] _, Multi-parameter paraproducts. Rev. Mat. Iberoam. 22 (2006), no. 3, 

963-976. T35 

[26] Zeev Nehari, On bounded bilinear forms, Ann. of Math. (2) 65 (1957), 153-162. 
12 

[27] Yumeng Ou, A T{b) theorem on product spaces, arXiv: 1305.1691, to appear 
Trans. Amer. Math. Soc., 1-45. |25 

[28] _, Multi-parameter singular integral operators and representation theo¬ 

rem, arXiv: 1410.8055, 1-28. t4, 25, 33, 34 


37 








[29] Stefanie Petermichl, Dyadic Shifts and a Logarithmic Estimate for Hankel Op¬ 
erators with Matrix Symbol, Comptes Rendus Acad. Sci. Paris 1 (2000), no. 1, 
455-460. 124 

[30] Stefanie Petermichl, Sergei Treil, and Alexander Volberg, Why the Riesz trans¬ 
forms are averages of the dyadic shifts?, Publicacions matematiques, ISSN 
0214-1493 46 (2002), no. 2, 209-228. 124 

[31] Akihito Uchiyama, A constructive proof of the Fefferman Stein decomposition 
of BMO{MA)., Acta Math. 148 (1982), 215-241. 13 


38 


